% Written by: Orhan Torul
% Istanbul, Turkey, May 2017
function Gvec = sguhm(varhm,fy1,fy,fx1,fx)

global ga ha;

a11_g=varhm(1);
a12_g=varhm(2);
a21_g=varhm(3);
a22_g=varhm(4);
a31_g=varhm(5);
a32_g=varhm(6);
a41_g=varhm(7);
a42_g=varhm(8);
a51_g=varhm(9);
a52_g=varhm(10);
a61_g=varhm(11);
a62_g=varhm(12);
a71_g=varhm(13);
a72_g=varhm(14);
a81_g=varhm(15);
a82_g=varhm(16);
a11_h=varhm(17);
a12_h=varhm(18);
a21_h=varhm(19);
a22_h=varhm(20);
f_y_1=fy1;
f_y=fy;
f_x_1=fx1;
f_x=fx;
ga=[a11_g a12_g; a21_g a22_g; a31_g a32_g; a41_g a42_g; a51_g a52_g; a61_g a62_g; a71_g a72_g; a81_g a82_g];
ha=[a11_h a12_h; a21_h a22_h];

Gara =f_y_1*ga*ha + f_y*ga + f_x_1*ha + f_x;
Gvec=Gara(:);
g_matrix_hm=ga;
h_matrix_hm=ha;
